Transposable element annotation in non-model species: The benefits of species-specific repeat libraries using semi-automated EDTA and DeepTE de novo pipelines.

Transposable elements (TEs) are significant genomic components which can be detected either through sequence homology against existing databases or de novo, with the latter potentially reducing the risk of underestimating TE abundance. Here, we describe the semi-automated generation of a de novo TE library using the newly developed EDTA pipeline and DeepTE classifier in a non-model teleost (Corydoras fulleri). Using both genomic and transcriptomic data, we assess this de novo pipeline's performance across four TE based metrics: (i) abundance, (ii) composition, (iii) fragmentation, and (iv) age distributions. We then compare the results to those found when using a curated teleost library (Danio rerio). We identify quantitative differences in these metrics and highlight how TE library choice can have major impacts on TE-based estimates in non-model species.

Ageing European lobsters (Homarus gammarus) using DNA methylation of evolutionarily conserved ribosomal DNA.

Crustaceans are notoriously difficult to age because of their indeterminate growth and the moulting of their exoskeleton throughout life. The poor knowledge of population age structure in crustaceans therefore hampers accurate assessment of population dynamics and consequently sustainable fisheries management. Quantification of DNA methylation of the evolutionarily conserved ribosomal DNA (rDNA) may allow for age prediction across diverse species. However, the rDNA epigenetic clock remains to be tested in crustaceans, despite its potential to inform both ecological and evolutionary understanding, as well as conservation and management practices. Here, patterns of rDNA methylation with age were measured across 5154 bp of rDNA corresponding to 355 quality-filtered loci in the economically important European lobster (Homarus gammarus). Across 0- to 51-month-old lobsters (n = 155), there was a significant linear relationship between age and percentage rDNA methylation in claw tissue at 60% of quality-filtered loci (n = 214). An Elastic Net regression model using 46 loci allowed for the accurate and precise age estimation of individuals (R 2 = 0.98; standard deviation = 1.6 months). Applying this ageing model to antennal DNA from wild lobsters of unknown age (n = 38) resulted in predicted ages that are concordant with estimates of minimum size at age in the wild (mean estimated age = 40.1 months; range 32.8-55.7 months). Overall, the rDNA epigenetic clock shows potential as a novel, nonlethal ageing technique for European lobsters. However, further validation is required across a wider range of known-age individuals and tissue types before the model can be used in fisheries management.

Contemporary evolution of the innate immune receptor gene TLR3 in an isolated vertebrate population.

Understanding where genetic variation exists, and how it influences fitness within populations is important from an evolutionary and conservation perspective. Signatures of past selection suggest that pathogen-mediated balancing selection is a key driver of immunogenetic variation, but studies tracking contemporary evolution are needed to help resolve the evolutionary forces and mechanism at play. Previous work in a bottlenecked population of Seychelles warblers (Acrocephalus sechellensis) show that functional variation has been maintained at the viral-sensing Toll-like receptor 3 (TLR3) gene, including one nonsynonymous SNP, resulting in two alleles. Here, we characterise evolution at this TLR3 locus over a 25-year period within the original remnant population of the Seychelles warbler, and in four other derived, populations. Results show a significant and consistent temporal decline in the frequency of the TLR3C allele in the original population, and that similar declines in the TLR3C allele frequency occurred in all the derived populations. Individuals (of both sexes) with the TLR3CC genotype had lower survival, and males - but not females - that carry the TLR3C allele had significantly lower lifetime reproductive success than those with only the TLR3A allele. These results indicate that positive selection on the TLR3A allele, caused by an as yet unknown agent, is driving TLR3 evolution in the Seychelles warbler. No evidence of heterozygote advantage was detected. However, whether the positive selection observed is part of a longer-term pattern of balancing selection (through fluctuating selection or rare-allele advantage) cannot be resolved without tracking the TLR3C allele over an extended time period.

Removal of beneficial insertion effects prevent the long-term persistence of transposable elements within simulated asexual populations.

BACKGROUND: Transposable elements are significant components of most organism's genomes, yet the reasons why their abundances vary significantly among species is poorly understood. A recent study has suggested that even in the absence of traditional molecular evolutionary explanations, transposon proliferation may occur through a process known as 'transposon engineering'. However, their model used a fixed beneficial transposon insertion frequency of 20%, which we believe to be unrealistically high. RESULTS: Reducing this beneficial insertion frequency, while keeping all other parameters identical, prevented transposon proliferation. CONCLUSIONS: We conclude that the author's original findings are better explained through the action of positive selection rather than 'transposon engineering', with beneficial insertion effects remaining important during transposon proliferation events.

No severe genetic bottleneck in a rapidly range-expanding bumblebee pollinator.

Genetic bottlenecks can limit the success of populations colonizing new ranges. However, successful colonizations can occur despite bottlenecks, a phenomenon known as the genetic paradox of invasion. Eusocial Hymenoptera such as bumblebees (Bombus spp.) should be particularly vulnerable to genetic bottlenecks, since homozygosity at the sex-determining locus leads to costly diploid male production (DMP). The Tree Bumblebee (Bombus hypnorum) has rapidly colonized the UK since 2001 and has been highlighted as exemplifying the genetic paradox of invasion. Using microsatellite genotyping, combined with the first genetic estimates of DMP in UK B. hypnorum, we tested two alternative genetic hypotheses ('bottleneck' and 'gene flow' hypotheses) for B. hypnorum's colonization of the UK. We found that the UK population has not undergone a recent severe genetic bottleneck and exhibits levels of genetic diversity falling between those of widespread and range-restricted Bombus species. Diploid males occurred in 15.4% of reared colonies, leading to an estimate of 21.5 alleles at the sex-determining locus. Overall, the findings show that this population is not bottlenecked, instead suggesting that it is experiencing continued gene flow from the continental European source population with only moderate loss of genetic diversity, and does not exemplify the genetic paradox of invasion.

Satyrization in Drosophila fruitflies.

The satyr of Greek mythology was half-man, half-goat, with an animal persona signifying immoderate sexual appetites. In biology, satyrization is the disruption of reproduction in matings between closely related species. Interestingly, its effects are often reciprocally asymmetric, manifesting more strongly in one direction of heterospecific mating than the other. Heterospecific matings are well known to result in female fitness costs due to the production of sterile or inviable hybrid offspring and can also occur due to reduced female sexual receptivity, lowering the likelihood of any subsequent conspecific matings. Here we investigated the costs and mechanisms of satyrization in the Drosophila melanogaster species subgroup of fruitflies. The results showed that D. simulans females experienced higher fitness costs from a loss of remating opportunities due to significantly reduced post-mating sexual receptivity than did D. melanogaster females, as a result of reciprocal heterospecific matings. Reciprocal tests of the effects of male reproductive accessory gland protein (Acp) injections on female receptivity in pairwise comparisons between D. melanogaster and five other species within the melanogaster species subgroup revealed significant post-mating receptivity asymmetries. This was due to variation in the effects of heterospecific Acps within species with which D. melanogaster can mate, and significant but nonasymmetric Acp effects in species with which it cannot. We conclude that asymmetric satyrization due to post-mating effects of Acps may be common among diverging and hybridising species. The findings are of interest in understanding the evolution of reproductive isolation and species divergence.

Epistatic regulation of growth in Atlantic salmon revealed: a QTL study performed on the domesticated-wild interface.

BACKGROUND: Quantitative traits are typically considered to be under additive genetic control. Although there are indications that non-additive factors have the potential to contribute to trait variation, experimental demonstration remains scarce. Here, we investigated the genetic basis of growth in Atlantic salmon by exploiting the high level of genetic diversity and trait expression among domesticated, hybrid and wild populations. RESULTS: After rearing fish in common-garden experiments under aquaculture conditions, we performed a variance component analysis in four mapping populations totaling ~ 7000 individuals from six wild, two domesticated and three F1 wild/domesticated hybrid strains. Across the four independent datasets, genome-wide significant quantitative trait loci (QTLs) associated with weight and length were detected on a total of 18 chromosomes, reflecting the polygenic nature of growth. Significant QTLs correlated with both length and weight were detected on chromosomes 2, 6 and 9 in multiple datasets. Significantly, epistatic QTLs were detected in all datasets. DISCUSSION: The observed interactions demonstrated that the phenotypic effect of inheriting an allele deviated between half-sib families. Gene-by-gene interactions were also suggested, where the combined effect of two loci resulted in a genetic effect upon phenotypic variance, while no genetic effect was detected when the two loci were considered separately. To our knowledge, this is the first documentation of epistasis in a quantitative trait in Atlantic salmon. These novel results are of relevance for breeding programs, and for predicting the evolutionary consequences of domestication-introgression in wild populations.

Comparative genetic stock structure in three species of commercially exploited Indo-Malay Carangidae (Teleosteii, Perciformes).

We examine genetic structuring in three commercially important species of the teleost family Carangidae from Malaysian waters: yellowtail scad Atule mate, bigeye scad Selar crumenophthalmus and yellowstripe scad Selaroides leptolepis, from the Indo-Malay Archipelago. In view of their distribution across contrasting habitats, we tested the hypothesis that pelagic species display less genetic divergence compared with demersal species, due to their potential to undertake long-distance migrations in oceanic waters. To evaluate population genetic structure, we sequenced two mitochondrial (mt)DNA [650 bp of cytochrome oxidase I (coI), 450 bp of control region (CR)] and one nuclear gene (910 bp of rag1) in each species. One hundred and eighty samples from four geographical regions within the Indo-Malay Archipelago including a population of yellowtail from Kuwait were examined. Findings revealed that the extent of genetic structuring among populations in the semi-pelagic and pelagic, yellowtail and bigeye were lower than demersal yellowstripe, consistent with the hypothesis that pelagic species display less genetic divergence compared with demersal species. The yellowtail phylogeny identified three distinct clades with bootstrap values of 86%-99% in mtDNA and 63%-67% in rag1. However, in bigeye, three clades were also observed from mtDNA data while only one clade was identified in rag1 dataset. In yellowstripe, the mtDNA tree was split into three closely related clades and two clades in rag1 tree with bootstraps value of 73%-99% and 56% respectively. However, no geographic structure appears in both mtDNA and rag1 datasets. Hierarchical molecular variance analysis (AMOVA), pair wise FST comparisons and the nearest-neighbour statistic (Snn ) showed significant genetic differences among Kuwait and Indo-Malay yellowtail. Within the Indo-Malay Archipelago itself, two distinct mitochondrial lineages were detected in yellowtail suggesting potential cryptic species. Findings suggests varying degrees of genetic structuring, key information relevant to management of exploited stocks, though more rapidly evolving genetic markers should be used in future to better delimit the nature and dynamics of putative stock boundaries.

Help or hindrance? The evolutionary impact of whole-genome duplication on immunogenetic diversity and parasite load.

Whole-genome duplication (WGD) events occur in all kingdoms and have been hypothesized to promote adaptability. WGDs identified in the early history of vertebrates, teleosts, and angiosperms have been linked to the large-scale diversification of these lineages. However, the mechanics and full outcomes of WGD regarding potential evolutionary impacts remain a topic of debate. The Corydoradinae are a diverse subfamily of Neotropical catfishes with over 170 species described and a history of WGDs. They are divided into nine mtDNA lineages, with species coexisting in sympatric-and often mimetic-communities containing representatives of two or more of the nine lineages. Given their similar life histories, coexisting species of Corydoras might be exposed to similar parasite loads and because of their different histories of WGD and genome size they provide a powerful system for investigating the impacts of WGD on immune diversity and function in an animal system. Here, we compared parasite counts and the diversity of the immune-related toll-like receptors (TLR) in two coexisting species of Corydoras catfish (C. maculifer and C. araguaiaensis), one diploid and one putative tetraploid. In the putative tetraploid C. araguaiaensis, we found significantly lower numbers of parasites and significantly higher diversity (measured by both synonymous and nonsynonymous SNP counts) in two TLR genes than in the diploid C. maculifer. These results provide insight into how WGD may impact evolution, in this case by providing greater immunogenetic diversity.

Evolutionary drivers of kype size in Atlantic salmon (Salmo salar): domestication, age and genetics.

The diversity of reproduction and associated mating patterns in Atlantic salmon (Salmo salar) has long captivated evolutionary biologists. Salmo salar exhibit strategies involving migration, bold mating behaviours and radical morphological and physiological change. One such radical change is the elongation and curvature of the lower jaw in sexually mature males into a hook-like appendage called the kype. The kype is a secondary sexual characteristic used in mating hierarchies and a prime candidate for sexual selection. As one of the core global aquaculture fish species, however, mate choice, and thus sexual selection, has been replaced by industrial artificial fertilization seeking to develop more commercially viable strains. Removal of mate choice provides a unique opportunity to examine the kype over successive generations in the absence of sexual selection. Here we use a large-scale common-garden experiment, incorporating six experimental strains (wild, farmed and wild × farmed hybrids), experiencing one to three sea winters, to assess the impact of age and genetic background. After controlling for allometry, fork length-adjusted kype height (AKH) was significantly reduced in the domesticated strain in comparison to two wild strains. Furthermore, genetic variation at a locus on linkage group SSA1 was associated with kype height, and a locus on linkage group SSA23 was associated with fork length-adjusted kype length (AKL). The reduction in fork length-AKH in domesticated salmon suggests that the kype is of importance in mate choice and that it has decreased due to relaxation of sexual selection. Fork length-AKL showed an increase in domesticated individuals, highlighting that it may not be an important cue in mate choice. These results give us insight into the evolutionary significance of the kype, as well as implications of genetic induced phenotypic change caused by domesticated individuals escaping into the natural environment.

Cryptic diets of forage fish: jellyfish consumption observed in the Celtic Sea and western English Channel.

To establish if fishes' consumption of jellyfish changes through the year, we conducted a molecular gut-content assessment on opportunistically sampled species from the Celtic Sea in October and compared these with samples previously collected in February and March from the Irish Sea. Mackerel Scomber scombrus were found to feed on hydrozoan jellyfish relatively frequently in autumn, with rare consumption also detected in sardine Sardina pilchardus and sprat Sprattus sprattus. By October, moon jellyfish Aurelia aurita appeared to have escaped predation, potentially through somatic growth and the development of stinging tentacles. This is in contrast with sampling in February and March where A. aurita ephyrae were heavily preyed upon. No significant change in predation rate was observed in S. sprattus, but jellyfish predation by S. scombrus feeding in autumn was significantly higher than that seen during winter. This increase in consumption appears to be driven by the consumption of different, smaller jellyfish species than were targeted during the winter.

Where are they all from? - sources and sustainability in the ornamental freshwater fish trade.

The global trade in ornamental fish involves c. 125 countries worldwide and is worth c. US $15-30 billion each year. This total is dominated (90%) by freshwater fishes, most of which are sourced from breeding facilities located in developing countries, typically in Asia or South America, but also in Israel, USA and Europe. Some fish are still obtained from natural (wild) sources, but the exact percentage of wild-caught fish is difficult to quantify given a lack of reliable data. Although c. 1000 species of freshwater fishes are widely available (from a total of > 5300 on sale), the most dominant freshwater fishes in the market comprise only 30 species from the orders Cyprinodontiformes, Perciformes, Characiformes and Siluriformes. In this perspectives review paper, illustrative example case studies of wild-fish collecting (Barcelos and Rio Xingu, Brazil) and breeding projects (Java, Indonesia) are described. In addition, wild-collecting expeditions to West Papua, Indonesia are discussed, focused on discovering novel species of rainbowfish (Melanotaeniidae) for breeding in captivity. Sustainability of the aquarium industry is considered in its broadest sense. The aquarium industry has been portrayed as both a major threat to natural ecosystems, but also as being part of the solution in terms of helping to maintain species when they have gone extinct in the wild or offering an income to impoverished citizens who might otherwise engage in much more destructive practices.

OSF-Builder: A New Tool for Constructing and Representing Evolutionary Histories Involving Introgression.

Introgression is an evolutionary process which provides an important source of innovation for evolution. Although various methods have been used to detect introgression, very few methods are currently available for constructing evolutionary histories involving introgression. In this article, we propose a new method for constructing such evolutionary histories whose starting point is a species forest (consisting of a collection of lineage trees, usually arising as a collection of clades or monophyletic groups in a species tree), and a gene tree for a specific allele of interest, or allele tree for short. Our method is based on representing introgression in terms of a certain "overlay" of the allele tree over the lineage trees, called an overlaid species forest (OSF). OSFs are similar to phylogenetic networks although a key difference is that they typically have multiple roots because each monophyletic group in the species tree has a different point of origin. Employing a new model for introgression, we derive an efficient algorithm for building OSFs called OSF-Builder that is guaranteed to return an optimal OSF in the sense that the number of potential introgression events is minimized. As well as using simulations to assess the performance of OSF-Builder, we illustrate its use on a butterfly data set in which introgression has been previously inferred. The OSF-Builder software is available for download from https://www.uea.ac.uk/computing/software/OSF-Builder.

How quantitative is metabarcoding: A meta-analytical approach.

Metabarcoding has been used in a range of ecological applications such as taxonomic assignment, dietary analysis and the analysis of environmental DNA. However, after a decade of use in these applications there is little consensus on the extent to which proportions of reads generated corresponds to the original proportions of species in a community. To quantify our current understanding, we conducted a structured review and meta-analysis. The analysis suggests that a weak quantitative relationship may exist between the biomass and sequences produced (slope = 0.52 ± 0.34, p < 0.01), albeit with a large degree of uncertainty. None of the tested moderators, sequencing platform type, the number of species used in a trial or the source of DNA, were able to explain the variance. Our current understanding of the factors affecting the quantitative performance of metabarcoding is still limited: additional research is required before metabarcoding can be confidently utilized for quantitative applications. Until then, we advocate the inclusion of mock communities when metabarcoding as this facilitates direct assessment of the quantitative ability of any given study.

Neutral variation does not predict immunogenetic variation in the European grayling (Thymallus thymallus)-implications for management.

Preservation of genetic diversity is critical to successful conservation, and there is increasing demand for the inclusion of ecologically meaningful genetic information in management decisions. Supportive breeding programmes are increasingly implemented to combat declines in many species, yet their effect on adaptive genetic variation is understudied. This is despite the fact that supportive breeding may interfere with natural evolutionary processes. Here, we assessed the performance of neutral and adaptive markers (major histocompatibility complex; MHC) to inform management of European grayling (Thymallus thymallus), which routinely involves supplementation of natural populations with hatchery-reared fish (stocking). This study is the first to characterize MH II DAA and DAB loci in grayling and to investigate immune genetic variation in relation to management practice in this species. High-throughput Illumina sequencing of "introduced," "stocked native" and "non-stocked native" populations revealed significantly higher levels of allelic richness and heterozygosity for MH markers than microsatellites exclusively in non-stocked native populations. Likewise, significantly lower differentiation at the MH II than for microsatellites was apparent when considering non-stocked native populations, but not stocked populations. We developed a simulation model to test the effects of relaxation of selection during the early life stage within captivity. Dependent on the census population size and stocking intensity, there may be long-term effects of stocking on MH II, but not neutral genetic diversity. This is consistent with our empirical results. This study highlights the necessity for considering adaptive genetic variation in conservation decisions and raises concerns about the efficiency of stocking as a management practice.

Whole genome duplication and transposable element proliferation drive genome expansion in Corydoradinae catfishes.

Genome size varies significantly across eukaryotic taxa and the largest changes are typically driven by macro-mutations such as whole genome duplications (WGDs) and proliferation of repetitive elements. These two processes may affect the evolutionary potential of lineages by increasing genetic variation and changing gene expression. Here, we elucidate the evolutionary history and mechanisms underpinning genome size variation in a species-rich group of Neotropical catfishes (Corydoradinae) with extreme variation in genome size-0.6 to 4.4 pg per haploid cell. First, genome size was quantified in 65 species and mapped onto a novel fossil-calibrated phylogeny. Two evolutionary shifts in genome size were identified across the tree-the first between 43 and 49 Ma (95% highest posterior density (HPD) 36.2-68.1 Ma) and the second at approximately 19 Ma (95% HPD 15.3-30.14 Ma). Second, restriction-site-associated DNA (RAD) sequencing was used to identify potential WGD events and quantify transposable element (TE) abundance in different lineages. Evidence of two lineage-scale WGDs was identified across the phylogeny, the first event occurring between 54 and 66 Ma (95% HPD 42.56-99.5 Ma) and the second at 20-30 Ma (95% HPD 15.3-45 Ma) based on haplotype numbers per contig and between 35 and 44 Ma (95% HPD 30.29-64.51 Ma) and 20-30 Ma (95% HPD 15.3-45 Ma) based on SNP read ratios. TE abundance increased considerably in parallel with genome size, with a single TE-family (TC1-IS630-Pogo) showing several increases across the Corydoradinae, with the most recent at 20-30 Ma (95% HPD 15.3-45 Ma) and an older event at 35-44 Ma (95% HPD 30.29-64.51 Ma). We identified signals congruent with two WGD duplication events, as well as an increase in TE abundance across different lineages, making the Corydoradinae an excellent model system to study the effects of WGD and TEs on genome and organismal evolution.

Evolution: Fangtastic Venoms Underpin Parasitic Mimicry.

Venomous teeth are rare in fishes, which typically utilise spines for defence. A new study reveals the evolutionary origins of fangs and venom in the Nemophini blennies and shows that, in contrast to snakes and lizards, the fangs pre-date the venom.

Jellyfish on the menu: mtDNA assay reveals scyphozoan predation in the Irish Sea.

Localized outbreaks of jellyfish, known as blooms, cause a variety of adverse ecological and economic effects. However, fundamental aspects of their ecology remain unknown. Notably, there is scant information on the role jellyfish occupy in food webs: in many ecosystems, few or no predators are known. To identify jellyfish consumers in the Irish Sea, we conducted a molecular gut content assessment of 50 potential predators using cnidarian-specific mtDNA primers and sequencing. We show that jellyfish predation may be more common than previously acknowledged: uncovering many previously unknown jellyfish predators. A substantial proportion of herring and whiting were found to have consumed jellyfish. Rare ingestion was also detected in a variety of other species. Given the phenology of jellyfish in the region, we suggest that the predation was probably targeting juvenile stages of the jellyfish life cycle.

Does density influence relative growth performance of farm, wild and F1 hybrid Atlantic salmon in semi-natural and hatchery common garden conditions?

The conditions encountered by Atlantic salmon, Salmo salar L., in aquaculture are markedly different from the natural environment. Typically, farmed salmon experience much higher densities than wild individuals, and may therefore have adapted to living in high densities. Previous studies have demonstrated that farmed salmon typically outgrow wild salmon by large ratios in the hatchery, but these differences are much less pronounced in the wild. Such divergence in growth may be explained partly by the offspring of wild salmon experiencing higher stress and thus lower growth when compared under high-density farming conditions. Here, growth of farmed, wild and F1 hybrid salmon was studied at contrasting densities within a hatchery and semi-natural environment. Farmed salmon significantly outgrew hybrid and wild salmon in all treatments. Importantly, however, the reaction norms were similar across treatments for all groups. Thus, this study was unable to find evidence that the offspring of farmed salmon have adapted more readily to higher fish densities than wild salmon as a result of domestication. It is suggested that the substantially higher growth rate of farmed salmon observed in the hatchery compared with wild individuals may not solely be caused by differences in their ability to grow in high-density hatchery scenarios.